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Abstract: We consider interjet observables in hard QCD processes given by the 
energy flow E out in a region away from all hard jets. Here the QCD radiation is 
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non-perturbative effects. We derive an evolution equation (in the large N c limit) 
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and identify universal and geometry-dependent contributions. Our analysis confirms 
features found using numerical methods by Dasgupta and Salam. 
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1. Introduction 

Study of "interjet" hadronic emission [1, 2, 3] is of special interest in QCD. In general, 
such emission originates from the flow of colour between jets, and therefore its analy- 
sis is important for understanding the mechanism of the overall colour neutralization. 
In practical terms, the interjet radiation is typically soft and so the distributions are 
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sensitive to possible non-perturbative components such as underlying events. For 
recent and not-so- recent experimental studies in hadronic colliders, see [4, 5]. 

A characteristic feature of inclusive interjet distributions is that collinear singu- 
larities do not play a relevant role and that the distributions are essentially deter- 
mined by large angle soft emissions, which are responsible for the coherence of QCD 
radiation [6, 7]. 

An example of an interjet quantity is E ont , the total energy (or transverse mo- 
mentum) of hadrons emitted in a region C out away from all hard jets. The E out 
distribution is infrared (and collinear) safe so that all its perturbative (PT) coef- 
ficients are finite and computable (in principle, although the resulting series then 
does not converge). Typically one has E ont <C Q, with Q the hard scale of the 
process, so that reliable QCD estimates require the resummation 1 of the logarith- 
mically enhanced terms generated by soft emitted or virtual partons. Since for this 
observable no logarithms are generated from collinear singularities, one has that the 
leading contributions to the distribution are given by single logarithmic (SL) terms 
<ln n Q/£ out . 

This should be contrasted with the case of observables in which the measured 
region includes one or more jets. An example of this is the thrust in the e + e~ process. 
In such cases both soft and collinear singularities contribute, so that the logarithm 
of the distribution contains, besides SL terms, also double logarithmic (DL) terms 
(a"L n+1 with L a large logarithm). 

Typically one expects that less singular distributions would be more difficult to 
resum. This is actually the case here. For a distribution in a "global" observable, i.e. 
involving the entire phase space, the PT expansion can be resummed [8] by means 
of a linear evolution equation (a generalization of the DGLAP equation [9]) based 
on the factorization of collinear singularities and coherence of the QCD branching 
structure [6, 7]. The resulting distribution can be expressed as a Sudakov form 
factor (the exponential of a radiator) which corresponds to bremsstrahlung emission 
directly from the primary partons. 

As shown by Dasgupta and Salam [10], "non-global" observables 2 , i.e. involving 
a part of the phase space, cannot be described only by the bremsstrahlung process 
but they involve the entire structure of successive parton branching. Calculations in 
e + e~ and DIS were performed numerically by a Monte Carlo method based on dipole 
branching emission derived from the distribution of many soft gluons emitted off a 
colour singlet dipole of hard partons given in [7, 11] in the large N c limit. 

Given the relevance of interjet observables, it would be important to have an 
analytic formulation for the SL resummation and this is the aim of the present study. 

1 Only the average value of E ou t or the distribution for relatively large E ut are obtained by 
finite order calculations [1] . 

2 Interjet observables are a particular case of non-global observable. Their numerical study for 
e + e" has been performed in [2]. 
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We derive an evolution equation based on the soft multi-gluon distribution given in 
[7, 11] valid for soft emissions in all angular regions. This evolution equation resums 
all SL contributions for the interjet distribution. It is based on energy ordering 
but implies also angular ordering [6, 7], which is the basis of Monte Carlo QCD 
simulations [12]. As in [2] this formulation of the branching involves colour singlet 
dipoles in the large N c approximation. 

In the main text we consider the case of e + e~ annihilation with an unobserved 
region C in defined by a cone around the thrust axis. Our analysis is extended in 
Appendix A to processes with incoming hadrons such as DIS and hadron-hadron 
collisions with large Pt-jets. We show that for all these hard processes the interjet 
distributions are given (in the large N c limit) in terms of a single function which 
depends on the geometry of the process. 

We derive the shape of the interjet distribution for large Q/E out and the structure 
of the branching and we discuss emission into the C ont region. We characterize these 
aspects in terms of quantities which are either universal or geometry-dependent. We 
confirm features found in [2] in numerical studies. 

The paper is organized as follows. In section 2 we describe the observable we 
consider and in section 3 we express the interjet distribution in terms of soft multi- 
parton emission. In section 4 we derive the evolution equation and we compare the 
case of the interjet distribution with the ones of other observables. In section 5 we 
describe the general features of the distribution, with both analytical and numerical 
methods. Section 6 contains the detailed technical discussion for large Q/E out . Uni- 
versal and geometry-dependent features are derived. Finally, in section 7 we recall 
the main physical features of the interjet distribution. 

2. The observable 

For a given hard process we consider as our interjet observable the total energy of 
hadrons emitted in a phase space region C out away from all hard jets 



Most of the energy flows inside the jet regions so that typically we have E ont <C Q with 
Q the hard scale. The features we will describe for this observable can be extended 
to other interjet observables such as the total hadron transverse energy. The interjet 
region C out could be defined in various ways according to the hard process. The 
complementary region including the jets will be denoted by C in . 
In e + e~ annihilation we will consider the interjet region C out as 




(2.1) 



h e c out 



— cos 9 hl < cos 9 h < cos 9 V 



(2.2) 
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with $h the angle with respect to the thrust axis Ht- We will discuss the following 
collinear and infrared safe distribution 




with da n the n-hadron distribution and &t the total cross section. The dependence 
on 6 m and Q is understood. This distribution is normalized to one at the kinematical 
limit E out ~ Q. The process is dominated by two-jet events aligned along the thrust 
axis so that E out is typically much smaller then Q. Three-jet events, with E out ~ Q, 
are of order of a s (Q). 

In Appendix A we discuss examples of a similar interjet observable in DIS and 
hadron-hadron collisions with high P r jets. 



3. QCD resummation 

At parton level, for small E out the distribution in e + e~ is described by the emission 
of the primary quark-antiquark pair pp accompanied by secondary soft gluons ki, 

e + e~ — > p ki . . . k n p . (3.1) 

In the soft limit, the primary quark and antiquark are aligned along the thrust axis. 
The leading contribution of the soft multi-parton distribution M% is obtained by 
considering strong energy ordering for the emitted partons. Here, the phase space 
integration can be approximated by (E ~ Q) 

d 2 Q 

d$ n = JJc^—ie^ - Ui) . (3.2) 

i 

For any one of the n\ strongly energy-ordered regions, the real emission contribution 
to the soft multi-parton distribution is given, in the large N c limit, by the factorized 
expression (see [7, 11]) 



M^ph . . . k n p) = M^(pp) ■ S pP (h ...k n ), 

Sppih . . . kn) = — Y[ "I ^2 W n{P k h ■ ■ ■ ki n P) , «s = N c — , 

I 1 7T n 



(3.3) 



where the sum is over all n\ permutations. For the fundamental permutation we have 

W ^ k ' ■ ■ ■ « = (pMfetl.fer) • (qq,) s 1 " c ° s "« ■ (3A) 

This distribution has both soft (ui — > 0) and collinear (6 qq > — > 0) singularities. It is 
valid in any strongly energy-ordered region. No collinear approximations are involved 
in the derivation, so it is valid even at large angles as needed for our study. 
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To leading order for small E out , the distribution £ e + e - (E out ) is obtained by using 

— = S pP {ki ...k n )- d® n . (3.5) 

Here we have to add to the soft distribution given in (3.3) the virtual corrections 
to the same order in the soft limit. This will be done in the next section in which 
we derive the evolution equation giving S e + e -(£ , out ). The expression (3.5), being 
accurate at the leading order in the soft limit, contains all the SL terms we want to 
resum. 

The basis for the resummation is the factorized structure of M%. We then fac- 
torize 3 also the theta function in (2.3) by writing 

\E ont - £>i = / 9 ■ I]*) ' < k ) = ©inW + e- w e out (A;) , 



iec, 



out 



2ttw 

(3.6) 

where O ut(&) and Oi n (&) are the support functions in the interjet region C out and 
Ci n respectively. The geometry of the interjet region C out , i.e. the 6 in dependence, is 
totally contained in the source u(k). The distribution is then given by 

/dv e vE ° ut _j 
— ^— — Gpp(E, v ) ~ G PP (E, E out ) , (3.7) 

where the Mellin variable v runs along the imaginary axis to the right of the singu- 
larities of G p p(E, z/ -1 ). In (3.7) we have evaluated the Mellin integration by steepest 
descent to give v ~ -Eout- The contribution from real emission (3.3) is given by 

G P T\E,E ont ) = 1 + jr [f[{a s ^^u(k t )Q(E-u t ))w n (pk 1 ...k nP ). (3.8) 

n=l J i=l < Ui 47F ) 

Here we set v = E~^ t in the source u(k) and used the symmetry in the secondary 
gluons to select the fundamental permutation at the cost of the l/n\ factor. Virtual 
corrections to (3.8) will be introduced in the next section. 

In Appendix A we extend this analysis to hard processes with incoming hadrons 
and show that interjet distributions for all hard processes are described by the same 
function G PiPj (E, E out ) with and pj the direction of pairs of hard jets (incoming 
or outgoing) and on the interjet region C ou t- 



4. Evolution equation 

To formulate an evolution equation we need to introduce the distribution G a b for a 
general pair of primary partons p a pb, obtained from (3.8) by replacing pp by generic 

Alternatively, we can use strong energy ordering so that E ont is the energy of the hardest 
secondary soft gluon. The Mellin transform method simplifies the combinatorics. 
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dipole momenta p a Pb- This distribution depends on the direction of p a Pb with respect 
to the e + e~ thrust axis and on the geometry of the interjet region C out , i.e. on 9\ n in 
(2.2). To obtain the evolution equation we use 



Ed E \ W{p a k x . . . k nPb ) J] e(E-Ui) \ = E5(E-u e ) w ab (k e 



i=l ) E=l 



(4.1) 



W( Pa h ...k e )H e(E-cji) \ ■ \ W(k e . . . k nPb ) 11 Q(E-uJi 

i=e+i 



where 



w (k) = ( PaPb ^ = 1 - cos 6 ab 

ab (p a k)(kp b ) (1 - cos6U)(l - cos6> fe6 ) ' 

We then deduce the basic equation {y = E~^ t dependence is understood) 

Ed E G ab {E) = J ^* a s w ab (k) [u(k) G ak (E) ■ G kb (E) - G ab {E)\ . (4.3) 

Here we have added virtual corrections, the last term in the square bracket, to 
the same order as the real emission contributions, see [7]. This evolution equation 
corresponds to soft dipole emission with energy ordering. Since w ab {k) effectively 
constrains k into the angular region within the ab dipole, (4.3) also implies angular 
ordering (after azimuthal averaging). Large angle regions are correctly taken into 
account. 

It is convenient to write (4.3) in the form 
E8 E G ab (E) = -Ed E R^h (E) ■ G ab (E) 

+ 1 ^ Wa ^ k) u ^ k) [Gafc(E) ' Gkb ^ E) ~ Gab ^ E) ] ' (4 ' 4) 

with (E) the SL Sudakov radiator for the bremsstrahlung emission 

Rf b \E) = J E ^.J ^a s w ab (k)[l-u(k)] = A-r ab , (4.5) 

where A depends on E, E out and r ab on the geometry of the interjet region (2.2) 

A = f E - a s [1 - e-^°-] , r ab = [ ^ w ab (k) . (4.6) 
Jo u Jc out 4vr 

Here we have used [1 — u(k)] = [1 — e~ w / E ° ut ] ou t(^) which entails that, in the un- 
observed jet region C in , the infrared and collinear singularities of w a b(k) are fully 
cancelled between real and virtual contributions. 
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In the quantity A the argument of the running coupling a s is determined by 
the successive hard emission contributions and cannot be determined in the present 
analysis. Detailed analysis [13] shows that, in the physical scheme, it is given by the 
transverse momentum relative to the afr-dipole which here can be approximated by 
the energy since the contributions to this radiator come from large angle emission. 

For small E out we can evaluate A by the standard SL approximation 

1 - e-"*" ~ Q(lu - e^v- 1 ) ~ Q(lu - E out ) , A ~ / — a s . (4.7) 

The radiator is then the contribution from the virtual parton in the C out phase 
space with energy u > E ont . For fixed a s we have A ~ a s In E/E out . 

The evolution equation (4.4) does not generate any collinear logarithms, either 
from the bremsstrahlung factor r ab (since the region C out does not include p a or p b ) 
or from the integral term in the second line of (4.4). The latter is collinear regular 
due to cancellation between the branching G ak • G kb and the virtual — G ab terms. 
Indeed, for k collinear to p a one has G a k 1 and Gkb G ab so that the sum in the 
square brackets vanishes and regularizes the singularity of w ab (k). Since only soft 
singularities remain, the leading terms of G ab are SL contributions. 

To SL order, we can replace u(k) by Q- m (k) in the integral term of (4.4). To show 
this observe that the remaining piece e~ w / E ° nt Q ou t(k) of the source contributes with 
energy uo < E out (see (4.7)) while the softest gluon, which contributes to R^j , has a 
typically larger energy (u > E out ). We conclude that the soft secondary branching 
takes place in C in and only the final parton enters the interjet region C out and here it 
contributes only with the virtual correction with energy oo > E out . 

By replacing u(k) — > Qm(k) in the branching term of (4.4) we have that the 
distribution G ab depends on E and E out only through the function A. We then can 
replace G ab (E, E out ) — > G ab (A) and obtain, to SL accuracy, 

d A G ab (A) = -r ab G ab (A) + I ^±w ab (k)[G ak (A).G kb (A)-G ab (A)\, (4.8) 

with the initial condition G ab (A — 0) — l. The distribution G ab (A) can be factorized 
into two pieces 

G afe (A)= e -^ ) ( A )^ a5 (A). (4.9) 

The first is the Sudakov factor given by bremsstrahlung emission from the primary 
hard partons p a p b . The second factor is the result of successive soft secondary branch- 
ing which satisfies 

d A g ab (A)=j c ^w ab (k) [u^liA) g ak (A) • g kb (A) - g ab (A)] , (4.10) 
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where 



U^ k (A) = e -^(««(A) , (4 . n) 

is an effective source. Notice that (4.10) has the same structure as (4.3) with u 
replaced by U^. The factorization in (4.9) can then be iterated, see later. 

4.1 Comparison with global observables 

The evolution equation (4.4) can be generalized to any collinear and infrared safe 
observable. Consider for example an additive global observable 

f = 5>(*i), (4-12) 

i 

with v(ki) a function of the emitted hadron momentum. One deduces 4 again (4.4) 
with the bremsstrahlung radiator involving the source for the specific observable 

R%\E) = J ^a sWab (k)[l-e-^]. (4.13) 

This quantity is collinear and infrared safe if v(ki) vanishes for ku — > 0. In this case 
the radiator involves also double logarithmic contributions a™ ln™ +1 V which arise 
from the collinear and infrared singular structure of w a b(k). SL terms coming from 
hard collinear emission are not included in (4.13), but these can easily be included 
by adding to w a f,(k) the finite pieces of the splitting functions. 

The soft secondary branching in the evolution equation (second line in (4.4)) 
plays a very different role 5 for global and local observables. As previously discussed, 
for local observables all soft secondary branching terms contribute to SL accuracy. 
For global observables, instead, soft secondary branching contributes beyond SL 
accuracy [6, 7]. This fact can be easily seen in the present formulation. 

Consider for instance the two loop contribution which involves the combination 
of the sources w(A; 2 )[l — u(ki)] with uj\ < cj 2 - For the local observable E out here 
considered, we have that w(/c 2 )[l — u(ki)] is consistent with energy ordering provided 
k\ G C out and k 2 G C in so that one has a SL contribution a 2 In 2 E/E ont . For a global 
observable instead, with C in = 0, we have that the combination -u(/c 2 )[l — u(ki)] 
contributes in the region v (k?) > v{k\) (see (4.7)) which is in conflict with energy or- 
dering. As a result, secondary branching has no infrared logarithms and one remains 
with a subleading term, a 2 In V, which comes from the collinear singularity of the 
dipole emission. As a result for a global observable, the distribution is simply given 
by a Sudakov form factor, the exponentiation of the DL bremsstrahlung radiator. 

4 In general one may need to factorize phase space constraints also, see [8]. 

5 As recalled before, the present treatment does not include hard contributions from the secondary 
branching which reconstruct the running coupling argument. 
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In (2.3) we have considered interjet observables in which C out does not include 
the primary parton region. The analysis can be generalized to the case in which the 
phase space region defining the local observable includes one primary parton [10]. 
Here one has DL contributions and soft secondary branching giving SL terms to all 
orders. 



5. General features of the distribution 

We cannot solve analytically the SL evolution equation (4.8). However we can study 
various aspects of the distribution using different approaches to give a picture of 
how the distribution behaves. We first observe that the distribution satisfies the 
constraint 



< G ab (A) < 1 , 



(5.1) 



for any A and ab. A proof of this is given in Appendix B. We also know that 
C a fe(0) = 1 and that, since the afr-dipole does not emit for a = b we have 



G ab (A) - 1 , 



for 



(5.2) 



for any A. For a ^ b, G ab (A) — > as A becomes large. This implies that the 
distribution G ab (A) has a peak at a = b, whose width decreases with increasing A. 
This peak governs the evolution of the entire distribution. 

We have used three approaches to study the distribution: an iterative solution 
valid at small A, the asymptotic behaviour at large A, and numerical solution. We 
discuss the results in the following subsections. 

5.1 Iterative solution at small A 

We observe that (4.10) for g ab is similar to (4.3) for G ab with u(k) replaced by Jjf b \ 
defined in (4.11). By iterating the procedure used to factorize the bremsstrahlung 
piece (see (4.9)) we obtain the general expansion 



ft*(i*)=exp -5>i?(A) 



71=1 



with the radiator components defined iteratively by 



R 



(n+l) 
ab 



U$ k (A) = e -<^ A )-^(A)+<; ; (A) 



(A) = 



dA' 



d 2 n k 



n-1 



JO 47T 

o L. ln 



™ ab (k) Uug(A') l-^(A') 



i=0 



(5.3) 



(5.4) 



We see immediately that at small A, R^(A) is of order A n+1 so that all terms 
contribute to SL accuracy. This is the crucial difference with the case for global 
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observables in which the soft secondary branching contributes beyond SL level, see 
discussion in subsection 4.1. 

In Appendix C we compute the first term for the e + e _ physical distribution 
R^(A). At small A we find (see [2]) 

fl M(A) = A^^-^M) +0 (A3), (5.5) 

while at large A it becomes 

B$(&) = AlnjAe^-^l - tan 4 %)) + O (1) . (5.6) 

As A increases, the higher contributions will also become important. However it 
is interesting to note the dependence on 6- m - even for quite large 8- m (of order 1) 
the distribution R^p is relatively insensitive to it, and in order to get a noticeable 
# in -dependence one has to move close to the limit 9 in = |. 

5.2 Asymptotic behaviour at large A 

The most physically interesting study is of the behaviour at large A, which is dis- 
cussed at length in section 6. Here we summarise the method and important results. 

To a first and quite crude approximation, we may treat G a b(A) at large A as 
being approximately 1 around a = b, and very small elsewhere, in other words having 
a peak of height 1 at a = b. This cuts off the integral in equation (4.8) for the 
secondary emission k near the emitters a or b. Physically this means that real gluons 
are only emitted close to one of the emitters; in the other regions only virtual gluons 
contribute. This is the same effect as described in [2] , where from numerical results it 
was seen that secondary gluons were radiated in a region around the primary partons, 
but that there also existed a large empty "buffer" region extending to the edge of 
C in . Importantly, this also means that the geometry of the C in region becomes less 
important, since what determines if a gluon is emitted or not in a particular direction 
is no longer whether that direction is in C in or C out , but if it is within the peak around 
one of the primary emitters or not. Thus we find that the leading large- A behaviour is 
independent of geometry. (Of course subleading terms are still geometry-dependent.) 

Of particular physical importance is the width of the distribution peak, which 
is a measure of the size of the buffer. We quantify this in terms of a critical angle 
6? cr i t (A) around parton a. Then we introduce a function h&(z) which for large A 
we take to give the shape of the peak with normalised width. We derive coupled 
evolution equations for # crit (A) and h&{z) which we evolve numerically to give what 
we assume to be stable limits for d\n9 clit /dA and h^z). At large A we then find 

# crit (A)~A a (C in )-e-i A , (5.7) 
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where c ~ 2.5 is a universal constant. The subleading behaviour, and in particular the 
normalization constant A a , is dependent on the geometry of C m and on the position 
of parton a, see section 6. Thus at large A the region within # crit (A) rapidly shrinks 
and the interjet region C in becomes more and more empty. This compares well with 
[2], who suggest from the numerical analysis the same functional form, with a value 
for c of about 3 (which is probably consistent within errors). 

Then using the asymptotic peak shape h^z) we obtain the large- A behaviour 
for the whole distribution. Consider the physical distribution for a general hard 
process given in terms of G PiPj (A) with pi and pj the hard primary jet directions 
which in general form an angle 0^. We derive the following Gaussian behaviour at 
large A 

Here c is the same universal constant in (5.7) and d a second universal constant. The 
functions and Aj, Xj depend on the geometry C m and the jet directions. They arise 
as integration constants in the large A approximation of the evolution equation. 

When the excluded jet region C in is a pair of small cones of angle 9 in centred on 
the jet directions PiPj the width of the peak becomes proportional to 9 in 

K(Cm)=X j (C in )=e in -X . (5.9) 

We obtain 

G ^ A K 2(i-i<y ) A ^ A >- »<A) = «- ,4 '(tfi«-*) 4 tf- <™> 

At large A the Gaussian behaviour coming from the soft secondary emission domi- 
nates over the bremsstrahlung contribution, which is the first factor. For small 9- ni 
the soft secondary branching distribution g(A) becomes fully # in -independent and 
direction-independent. The constants / , A arise one factor from each jet, since for 
small C; n the two jet contributions factorize, and by symmetry they are equal. Thus 
in any hard QCD process the distribution of energy emitted into the region excluding 
a narrow cone around each hard jet is determined for each colour flow by the usual 
bremsstrahlung radiator multiplied once for each emitting dipole by the universal 
Gaussian function g(A). 

If the weak dependence of R^p(A) on 6 m found in (5.6) is at all indicative of the 
general behaviour, we may hope that the small cone approximation (5.10) for g(A) 
could be valid even for quite large cones: this was indeed seen in [2]. This smooth 
behaviour is explained by the fact that for small 6- in the critical region 6* crit (A) inside 
which the branching develops scales proportionally to 9 in , see (5.9). 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
1-cos6 ab 



Figure 1: The distribution g a b for cos#; n = 
three different values of A. 

5.3 Numerical results 



i and 9 a = as a function of 1 — cos 6 a b for 



We try a numerical solution of (4.10) to give the quantitative behaviour of the sec- 
ondary branching distribution g> &(A) in e + e~ annihilation. We set a grid of 40 bins 
in and 80 bins in cos8 and evolve (4.10) numerically. Real- virtual cancellation is 
implemented by neglecting the contribution with k = a or k = b in (4.10). This may 
be a crude approximation, due to the sharp behaviour of the function g a b for a near 
b, and needs improvement in order to have complete control on the errors. For what 
concerns the A behaviour of the solution, we find substantial agreement with the 
results obtained in [2]. 

As stated in the previous subsection, it is crucial to understand the behaviour 
of the function g a b when we vary the opening angle between a and b. In figure 1 we 
plot g a b as a function of 1 — cos 6 a b for three different values of A. We choose to fix 
a along the thrust axis and let 6 a b run from to 6 m . This distribution starts from 1 
and then steeply decreases with increasing 6 a b- This behaviour confirms the presence 
of a peak for the g a b distribution near a = b which shrinks with increasing A. The 
form of this distribution for large A will be discussed in detail in the next section. 
From this figure we see that the asymptotic behaviour is already settled at A = 3 
where the distribution is negligible away from the peak. 

6. Study of the large-A behaviour 

We focus on e + e~ annihilation with C m given by a cone of size 8 m around the two 
jets along thrust axis, see (2.2). We then generalize this analysis to the distribution 
GpiPj required for other hard processes, in which the two jets are not back-to-back. 
To obtain G a b(A) in the large A limit we need first to study its behaviour near the 
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peak at a = b. Using this information we can obtain the asymptotic behaviour of G ab 
for any a, b. These two stages are described in the next two subsections. In a final 
subsection we consider the particular case of C in being a pair of small cones. 

6.1 Shape of the peak at a = b 

For each point a inside the jet region C in , we can measure the width of the peak in b 
near a by using 

f ^ G ab (A) = 1 - costf crit (A) « iCit(A) • (6.1) 

In general this quantity depends on the geometry of the C in region and on the point 
a chosen in C in . Here # crit (A) measures the angle between a and b above which the 
distribution G ab is suppressed. It is therefore a small angle which decreases as A 
increases. We have the initial value 6* cr i t (A = 0) = #; n . 

Since the evolution of the peak is determined only by the distribution around the 
peak, and not by points remote in phase-space, we have that the shape of the peak 
in G a {, (A) depends only on the angle 9 ab between a and b. To measure this shape, 
suppose that at large A and in the region near a = b the distribution behaves as 

G a6 (A) = h A (z) , z = , (6.2) 

for some function h A (z) that from (6.1) for all large A must satisfy 

dzh A {z) = l. (6.3) 



Jo 



The evolution equation (4.8) simplifies in that the source term involving r ab is neg- 
ligible and the integration over solid angle becomes an integration over a flat plane: 



(6.4) 



0aM*) " ' zh' A (z) = -L J [h A (rlz)h A (rjz) - h A (z)} , 

r\ = x 2 + y 2 , r\ = (1 - xf + y 2 , 9 ak = r x e ah , 9 kb = r 2 9 ab . 
Integrating this equation over z yields 

J_ f°° dz f [h A (r 2 lZ )h A (r 2 2 z) - h A {z)] . (6.5) 

^ J0 J r l r 2 



dln9 2 Tit 



dA 

Equations (6.3)-(6.5) form a coupled system that can be evolved numerically from 
any suitable starting function ho(z). Here h A (z) is a function of only one variable z 
while the full distribution G ab (A) is a function of three angles, so the implementation 
is somewhat easier. 
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There are many possible choices of starting function ho(z): it need only be 
bounded between and 1, with ho(0) = 1, normalised, and sufficiently smooth. We 
used a selection of such functions, namely 

ho(z) = e-H 1+ ^r ,(„=i, 1, 2, 3, 4) , h (z) = ^6(2 - z) , (6.6) 

all of which showed the same asymptotic behaviour. The evolution (6.4) settles down 
to a shape h^z) shown in figure 2, which satisfies the equation 



c - zh 'oo( z ) = 7T [ —TT [^00(^1^)^00(^2^) - hoo(*)] , 
Z7r j r 1 r 2 

~27r J dZ J [ h °°( r i z ) h °°( r 2 z ) ~ hoo( z )] ~ 2.5 , 



C= dA~~ 

with c evaluated with an accuracy of about 10%. This error is due to limitations in the 
numerical analysis from the finite number of ^-points, and the numerical integrals 
performed at each stage, such that evolving using the same starting function but 
with a different number of points or precision in the integrals converges to a slightly 
different numerical value for c. The step size in A needs to be small enough to prevent 
instabilities developing before convergence is seen — we found the value 0.02 to be 
sufficiently small. The stated uncertainty on the value of c is a generous estimate. 
This error could be improved with a more refined analysis. 

We then conclude that 9 CTit (A) has the behaviour at large A given in (5.7) with 
A a an integration constant which depends on 9 in and on the chosen point a. The 
fact that /loo i z ) is finite for any z implies that at large A the distribution G a b(A) for 
&ab 1 depends on A,0 in and the point a only through the function # crit (A). The 
function h^z) as well as the constant c is universal. 

From the tail of the function h oc (z) at large z we can estimate the large- A 
behaviour of G a b(A) away from the peak in the region 

crit (A) < 9 ab < 1 . (6.8) 

To obtain this we rewrite equation (6.7) in the form 
din hoo If dxdy 



din ' TT J ri<r2 r\r\ 



h C o(r 1 z) - 1 + h^ir^z) 



hoo(r 2 z) 



hoo(z) 

As z — > 00 the final term in the brackets vanishes and we obtain 

d In hoo 



(6.9) 



with 



dlnz 
1 

T 

fa) 



= f 7[M^')-ik(;) , (6.10) 



x < 7 

x 4 




2 tan- M±± , x> 1 



(6.11) 



n(l-x) \ 1 + y/x V 2^fx~- 1 / ' 4 
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Figure 2: The large A peak, h^z). 



In the large z limit the evolution equation becomes 

d In h C 1 dz' r°° dz' 

c-^p^^-lnz-c', d = / ^ [1 - hM)] ~ / ^M*0- (6-12) 
amz Jo z Ji z 

Therefore, using (5.7) and (6.2), we conclude that for a and b in the region (6.8) we 
have, in the large A limit, 

lnG a6 (A) * -^A 2 - f c' + In || ) A - (- + °- In || + 1 In 2 || V (6.13) 



2 ^a/ V c c 2A2 2c 2A^ 

The constants c' and c" can be determined by the function h OQ (z) shown in figure 2. 
This behaviour is valid provided a, b are away from the boundary of C m . 

In taking the large- ,2 limit, we have neglected terms that contribute at finite z 
in moving from (6.9) to (6.10) and from (6.10) to (6.12). The leading correction is 
0(z^ x ), and arises from the expansion of the £ function in (6.10). Thus the right 
hand side of (6.13) has a correction O (e _cA ). 

6.2 The distribution off the peak 

Using the fact that we know the form of G a b(A) near the peak at a = b, we now 
determine the large-A behaviour of G a b(A) for any a and b including the physical 
case in which a and b are along the two hard jet directions (6 a b ~ 1)- 

We find that the geometry dependence in the integration regions of the two terms 
in (4.8) cancels for large A. To obtain this result we write (4.8) in the form 

d A lnG ab = - j ^w ab (k) + ! ^w ab (k) _ . (6 .i4) 

JC out JC in \ ab / 
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Since the first term in the second integral is negligible for k not near a or b we may 
write 



a 1 r> f d 2 Vt k f d 2 Vt k {G„k-G 

d A lnG ab ~- / — — w ab {k) + / — — w ab {k) 



ak ' t - r fe6 j 



"ak>"kb^^ "ak 



a b 



where e is a small parameter (larger than # cr i t ) on which the final answer should 
not depend. Here we have neglected the ratio G a kGkb/G a b for 6 a k, Q b k ^ $crit which 
vanishes at large A. Due to the cancellation of C in and C out in the integration region, 
we have that, at large A, the distribution G ab (A) depends on the geometry and on 
the directions of a and b only through 9 crit a (A) and 9 crit &(A), the critical values in 
(5.7). 

We now show that the e-dependence in the various terms cancel. The first term 
of (6.15) gives the contribution with a cutoff around a and b: it is, for small e, 

, ^^).ln^^, (6.16) 

while the second and third terms give the contributions from k near a or b, which 
depend on the shape of the peak given above. We write 



/ 

Je 2 



f ^^/G^^Li f ^(G,-l)~-^-iln^ 
4< e 4* Wab[k) { G ab 1 )-2j 6l k {Gak 1} " 2 2%^ 



(6.17) 

with # cr it ia the critical angle for the peak around a and d the constant given in (6.12). 
Here corrections vanish in the large A limit. A similar equation holds for the integral 
around b. Thus the e dependence cancels and we obtain the evolution equation at 
large A: 

d A lnG ab ~ -cA - c' - In 1 ~ C °^ ab . (6.18) 
The asymptotic behaviour of G ab is therefore 

G ab (A) ~ e-§ A2 f AaA&e ' C/ ) A / a6 (0 in ) , (6.19) 

V 1 - COS d ah J 

where f ab is independent of A at large A. For 8 ab C 1 we have, see (6.13), 

7+^k + ^ 2 M)- (62o) 

Because of this structure, we have that secondary branching is almost collinear to 
the partons a, b initiating the branching: in fact no farther away than the width of 
the peak. This implies that as long as a, b are not close to the boundary of C in the 
geometry dependence enters only through the parameter A a in # crit . 
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6.3 Small C- in region 

Since the secondary branching is almost collinear to the hard primary partons it 
becomes interesting to study the limit 9 in — > 0. In particular we would like to see 
how the buffer region behaves when the jet region is squeezed. We consider the two 
cases of ab in the same and in opposite jet regions. 

Same jet region. We derive the evolution equation 



dAGab — — r ab ■ Gab + / T, — ~ — "TwF 7\T \Gak Gkb ~ G ab ] , 

'\k\<i ^ (a - k) 2 (k - b) 2 

b) 2 \ . e a _ j. ft b _ 

''"" = a*)(l-b 2 ) ' a= 6- na > b= 9- nb 




(6.21) 



using rescaled angular variables a, b. The contribution to the evolution of G a b from 
the branching in the opposite region vanishes for ft in — > 0. This is an aspect of 
coherence of the QCD radiation. We notice that the explicit in dependence has 
disappeared from the equation. This implies that # crit (A) scales with ft in so that 

Aa(C in ) = in • A(a) , (6.22) 

giving (5.9) in the case a = 0. We have then that the buffer region still expands by 
increasing A. By performing the same analysis as in subsection 6.1 we find, at large 
A, for ab in the same jet region but away from the peak, the behaviour 

W.^ ^^ pBi). (6.23) 

Opposite jet regions. We introduce now rescaled angular variables for the two 
back-to-back jets (left and right jet). For 6 m — > the radiator is 

r ah = -21n^ - iln(l-a 2 ) - iln(l-6 2 ) , d=^-n a , b = - — n b . 

II l fin fin 

(6.24) 

This explicit 9- m dependence implies that G a b(A) is given by 6* 2 n A (the bremsstrahlung 

— * 

contribution) times a function of the rescaled variables a, b, whereas the same jet 
distribution depends only on the rescaled quantities. 

The a and b contributions enter (6.24) independently. This is due to the fact 
that w a b(k) splits into a sum of right and left pieces. As a consequence the evolution 
in the right and left regions develops independently and G a b(A) factorizes 

/ft \ 2A 

G ah (A) ~ • F a R (A) ■ F 6 L (A) . (6.25) 
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Here F^- and F fc L satisfy the right and left evolutions and for the first we find 



9aF* = \ HI - a 2 ) • F a R + y if - F*] , (6.26) 

where is the distribution for a and A; in the same (right) jet region discussed 
above. Proceeding as in subsection 6.2 we find the following asymptotic behaviour 
for the distribution with ab in opposite jet regions 

G ab (A) ^ e-f- 2 ( ^MM^ j ' / (S) / ( g) , (6 . 27) 

with d the integration constant given in (6.12) and / an integration constant de- 
pending only on the rescaled variable. This expression shows that for small 9 in the 
function f a b(6m) in (6.19) factorizes and depends on the rescaled variables. The 
# in -dependence is the one given by the bremsstrahlung contribution. The secondary 
branching contribution does not depend on 9 in and is decreasing with A with a 
Gaussian behaviour. For the physical e + e~ distribution we set a 2 = b 2 = 0. 

A generalization of this result is to the case where the two jets piPj are not 
back-to-back but are set at an angle 0^. The asymptotic behaviour is given by the 
expression (6.27) with {2}~ A replaced by {1 — cos#jj}~ A and with 9 a the angle of 
a with respect to the jet pi and 9b the angle of b with p^. The physical distribution 
G PiPj (A) is again obtained by setting a 2 = b 2 = 0. 

7. Discussion 

We summarize here the essential physical points of our results on energy flow away 
from hard jets. The interjet distribution reduces to G p p(A) in (3.7) for e + e~, while 
for DIS and hadron-hadron collisions it is given in terms of G PiPj (A), see (A. 5) and 
(A. 11). These distributions depend on the geometry of C ont and on the jet directions. 

A characteristic of interjet observables is that the leading contribution to their 
distributions are SL, originating from soft emission at large angles. Collinear sin- 
gularities here are subleading. The contribution from bremsstrahlung emission from 
the primary hard partons does not give the full SL structure. Indeed we find (see 
(4.9)) 

G PiPj (A) = e- R( ^ ■ g PiPj (A) , (7.1) 

with i2piPj-(A) the SL Sudakov radiator for emission into C out . The additional SL fac- 
tor g PiPj is generated by soft secondary emission, which can be described by successive 
colour singlet soft dipole emission. This implies that, in the present formulation, we 
work in the large N c approximation. The soft secondary branching develops with 
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decreasing energy inside the unobserved jet region C m until a soft gluon enters C out . 
The emission is also effectively angular ordered. 

At large A the distribution G PiPj (A) decreases with a universal Gaussian be- 
haviour, see equation (5.8), with a process-independent coefficient c that we have 
estimated to be c « 2.5 with a 10% accuracy. Since the bremsstrahlung radiator 
is proportional to A we see that the large A behaviour is dominated by the soft 
secondary branching contribution g p . p .(A). 

The origin of the universal Gaussian behaviour in A is a consequence of the struc- 
ture of the branching. The secondary branching generating g PiPj (A) develops within 
a small cone # cr i t (A) around the hard primary partons pi and pj, which decreases ex- 
ponentially for large A according to (5.7), again governed by the universal coefficient 
c. It is the development of soft secondary branching in this peak region which gen- 
erates the Gaussian behaviour in (5.8). At large A the region within crit (A) shrinks 
generating an empty buffer region first observed in [2]. For the case of figure 1 we 
have estimated that the asymptotic regime sets in around A = 3. At this value the 
distribution is negligible away from the peak and approaches the asymptotic limit of 
figure 2. 

Taking the jet region C in to be an arrangement of small cones centred on each 
hard jet leads to a universal secondary emission function g(A) that depends neither 
on the size of the cones nor on the directions of the jets (see (5.10)). 

The distribution G PiPj (A) is a function of E,E out and the coupling a s through 
the quantity A defined in (4.6). This involves an integral over the running coupling 
at low scales, and so acquires non-perturbative power corrections. These corrections 
can be estimated by the same method [14] used for the standard shape variables in 
e + e~ and are expressed in term of a unique parameter A NP given by the integral over 
the running coupling in the low energy region. We obtain 

A ee l E ^a s (u) [1 - e^/Vt] ~ f d ^a s {u) + N c ^- . (7.2) 

JO & JE out V ^Amt 

The non-perturbative parameter A NP (for the normalization used here see [15]) is 
measured both at LEP [16] and at HERA [17]. 

In hadron-hadron collisions one expects [1, 5, 12] additional large contributions 
to interjet observables coming from a soft underlying event generated by the incoming 
hadrons. We can write 

E out = X>* + E%£ , (7-3) 

*e c out 

where the sum includes all partons produced in the hard process, and the quantity 
is the contribution of the soft underlying event. Assuming the underlying 
emission and hard emission factorize, the final answer will be given by the distribution 
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studied here with the replacement 



A(E, E out ) -> A(E, E out - E%£) . (7.4) 

Up to now in hadron-hadron collisions only the mean value of E ont has been investi- 
gated both theoretically and experimentally. One shows that the hard contribution, 
which can be reliably computed by fixed order QCD calculation 6 , does not fit the 
data [1, 4, 5] but requires a sizable E^ contribution. 

The resummed PT distribution here discussed provides a way to further study 
the underlying event, in particular its factorization properties and its size. 
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6 The PT calculation in [1] has been done only to af order, now it is possible to perform the 
calculation at order [18] so that one can control the QCD scales. 
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A. Interjet observables in DIS and hadron-hadron collisions 



We introduce first the physical observables and then the QCD resummation. 
A.l The observable 

An example of an interjet distribution in DIS is given by 



s »*^fc/Mi:/^ t«,- e 4 (a.i 



with da n / dx^dQ 2 the n-hadron distribution for fixed Bjorken variable x-q and Q 2 . 
The events are dominated by the incoming and the outgoing hard jet. In the Breit 
frame the situation is similar to e + e~ with both jets aligned along the beam direction. 
We may define the C out region by 

-riV < ^ tan | < ^ . (A.2) 
An example of an interjet distribution in hadron-hadron collisions with hard jets is 

*W*-> = {E / ^} _1 E/^ e( £out - E , (A,) 



out 



with da n /drjdP t the n-hadron distribution associated to the emission of a jet with 
rapidity rj and transverse momentum P t . The jet could be defined by a /c t -jet finding 
algorithm [19]. Here the event is dominated by the presence of an additional recoiling 
outgoing jet. The interjet region C out should be defined by avoiding the two incoming 
and the two outgoing hard jets. It is given for instance by fixing a region in the 
rapidity rjh and azimuthal angle (f>h of emitted hadrons with respect to the triggered 
jet, see [4, 5] . 

The QCD resummation discussed for e + e~ annihilation in section 3 can be gen- 
eralized to the above processes due to the factorization of the parton process (see for 
instance [20]) as we shall discuss. 

A.2 QCD resummation: DIS case 

To express the result for EDis(-Eout) we write the DIS cross section in the form 

da ~ f dx /^ 7 ^\ u ~ Q , (A.4) 



dx-sdQ 2 J dxidx-^dQ 2 

where da 1+pl ^ P2 {xi, u) is the factorized parton distribution including the parton 
density function at the hard scale u ~ Q and the squared matrix element of the hard 
vertex 7* +pi — * p 2 . Here X\ is the momentum fraction of the hard parton pi coming 
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into the hard vertex in the Breit frame. In the soft limit p± and p 2 correspond to the 
directions of the incoming and outgoing jets respectively. 

By performing the same analysis as done in section 3 we obtain 



J DIS 



The factorization scale remains at fi since no observation is made in the collinear 
region. Here G PlP2 (Q, E out ) is the distribution introduced in (3.8) with p\p 2 aligned, 
in the Breit frame, along the beam direction. 

A. 3 QCD resummation: hadron-hadron case 

The situation here is more involved since we have four jets (two outgoing and two 
incoming). There are various hard elementary processes 

P1+P2 -> P3 + Pi , (A. 6) 

according to the nature (q, q or g) of the four involved partons. They are: 

A q(p l ) + q'(p 2 ) -> q(p 3 ) + q'{pi) , 

B q(p l ) + q(p 2 ) -> q(p 3 ) + q(p A ) , ^ 
C q(p l ) + q(p 2 ) -> g(ps) +g(P4) , 

D ^(Pl) + 9(P2) -»• 5(Ps) + ^(P4) , 

together with those obtained by crossing transformations. 

Energy flow in interjet regions for this process has been studied in [3] as far as 
the SL bremsstrahlung radiator piece. To construct the SL contribution from the 
soft secondary branching we need to represent the soft n-gluon emission in each of 
the hard processes (A. 6) as a combination of colour singlet dipole emissions (3.3) 
from the four hard partons. The square amplitude for the emission in (A. 6) of a 
single soft gluon k has been written in [21]. Taking the large N c limit, this result can 
be expressed [22] as a sum of dipole emissions Wij(k) from pairs piPj of hard partons 
with the factorized structure 

Kmsu ,* - Km 3 u ■ E c flLu E ^ > E c m M = 1 > ( A - 8 ) 

where M^^^^ is the exact hard elementary 2^2 distribution for the process 
(A. 6) with fi the type of parton pj. The coefficients C/^/ 2 / 3 / 4 are functions of the 
kinematical invariants s— (P1+P2) 2 , t— (pi— P3) 2 , u— (pi—Pi) 2 and can be found 
in [22] . Here 7 represents the set of all colour-connected partons in the specific hard 
process and the last sum is over all colour connected pairs ij in the set 7. 
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The simplest process is qq' — > qq' with only (/1/4) and (f 2 f 3 ) colour connections. 
Here the coefficients are independent of the invariants and, from symmetry, we have 

E E w ^ k ) = w ^ k ) + ■ ( A - 9 ) 

7 fa") 67 

The most complex process is gg — > for which all possible colour connections 
contribute. 

The distribution for the emission in the hard process (A. 6) of n additional soft 
gluons in the strongly energy ordered region is obtained by performing the analysis 
done in [7, 11] in the large N c limit. Each additional soft gluon contributes to one of 
the hard dipole emissions as done in (3.3). This allows us to extend the analysis of 
section 3 and deduce the interjet E out distribution in hadron-hadron collisions with 
a hard jet of rapidity r\ and transverse momentum P t . 

Consider the hadronic cross section for this process which we write in the form 

(A10) 

where g^7i/ 2 -*/ 3 / 4 is the parton process integrand which includes parton density func- 
tions and hard parton distribution ^} 1 f 2 f^f A for the elementary process (A. 6). The 
normalized E out distribution in a given interjet region C out is then given by 

^fif 2 f 3 h( p t, Kut) = E ^1/2/3/4 II G PiPj( P t, E ut) , 

7 fa) e 7 

(A.ll) 

where G PtPj (P t , E out ) is the distribution introduced in (3.8) associated with the dipole 
Pi,Pj- In the soft limit we can take p\p 2 in the direction of the two incoming hadrons 
and p 3 p 4 in the direction of the two outgoing hard jets. We report the explicit form 
of 'Ef 1 f 2 f 3 f 4 (P t , E out ) for the elementary processes in (A. 7): 

H A Z A = h A (s,t,u) G l4 G 23 , 

H B S B = h B (s, i, u) G 14 G 23 + h B (s, u, t) G 13 G 24 , 

H c S c = h c (s, i, u) G 34 G 13 G 24 + h c (s, u, i) G 34 G 14 G 23 , (A. 12) 

H B S D = h B (s,t,u)G 12 G 24 G 43 G 31 + h D (s,u,t) G 12 G 23 G 34 G 41 
+ h B (u,t, s) G 14 G 42 G 23 G 31 , 
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h A (s,t,u) = g^ 
h B (s,t,u) = h A (s,t,u) + 2g*%H 

h c (s,t,u) = g^[^-^) 

h%s,t,u) = 2g* (£^)(l-^-f + i) 

Table 1: The functions /i( A > B >C,D) taken from Ref. [22]. 



with Gij = G PiPj (P t , E out ) and the functions /i( A > B > c ' D ) given in table 1. Here 

H A = h A {s,t,u), 

H B = h B (s,i,u)+h B (s,u,t) : 

H c = h c (s,i,u) +h c (s,u,t) , 

H B = h D (s,i,u) + h B (s,u,i) +h B (u,t,s) , 

are the square amplitudes for the 2 — > 2 elementary processes (A. 7). 

The structure in (A. 12) in terms of the factorized distributions is derived 
from eq. 66 of Ref. [22] in which we have identified 2Cp with Ca since we work in 
the large N c limit. Going beyond the large N c approximation requires the analysis of 
colour interference between different hard partons. At the moment this can be done 
only for the bremsstrahlung contribution, elegantly computed in [3], see also [23]. 

In conclusion, for all considered hard processes, the interjet distributions are 
described by the universal function G PiPj (E, E ont ) which depends on the two hard jet 
directions pi and pj (incoming or outgoing) and on the interjet region C ou t- 



B. Theorem: < G ab (A) < 1 

In this appendix we show that the solution of the differential equation (4.8) with 
the given boundary condition is bounded below and above by and 1. Thus our 
distribution is physically meaningful. 
Proof: 



1. G a fe(A) is a continuous function of A and of a and b. 

2. G ab (0) = lVa,6. 

3. If a = b then G ah (A) = 1 VA. 

4. Suppose 3 a', b' and A' with Gv b /(A') > 1. Then by continuity 3 a", b" and 
A" with < A" < A', G ab (A") < IV a, b, and a" ^ b", GW,(A") = 1, 
d A <w(A") > 0. 
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5. But by (4.8) we have <9a (A") < which gives a contradiction. Therefore 
G ab (A) < IV a, 6, A. 

6. Suppose now 3 a',b' and A' with G a i b >(A') < 0. Then by continuity 3 a",b" 
and A" with < A" < A', G ab (A") > 0Va,6, and a" ^ b", GW'(A") = 0, 
<9 A GW(A") < 0. 

7. But by (4.8) we have d/±G a » b "(A") > which gives a contradiction. Therefore 
G ab (A) >0Va,6,A. 



C. Iterative solution 



Here we calculate the first term RpJ (A) in the iterative solution described in section 
5.1. First we have that for any a, b the bremsstrahlung emission is given by 



s a cos 6 a + COS 9 

Tab = TT In Z 

2 cos 9 n — cos 9 



m s b cos 6 b + cos 6[ n s a + s b 

+ — In — H — x 

2 cos # b — cos 9 hl 4 



In 



( 1 + cos 2 # in ) (1 + cos 6 a cos 9 b ) — 2 cos # in (cos # a + cos 6 b ) — sin 2 9 in sin 9 a sin 6*;, COS (pab 
( 1 + cos 2 9- m ) ( 1 + cos 9 a cos 6*b) + 2 cos ^in (cos # a + cos 9 b ) — sin 2 # in sin 9 a sin # b cos 4> ab ' 

(C.l) 



where = 1 for cos6> fe > and Sk = —1 for cos9 k < 0. Therefore we obtain R^j in 
the following cases 

1 + cos 6> in 



i?g>(A) = A-ln 



1 — cos 9 m 



R (0) / _ A £fe j n (1 - COS gin) (COS gfc + cos 9 in ) 

pk [ ' ~ ' 2 11 (1 + cos in ) (cos 9 k - cos in ) 

1 + COS 6 , in )(cOS0fc + COS # in ) 



(C.2) 



d(«)/ A \ _ A £^ l n 

fcp_l ; 2 (1 -cos0 in )(cos0 fc -cos0 in ) 

Now using the definition (5.4) we explicitly calculate 

d 2 Q, 



i$(A) 



"PPM l-^(A') 



where 



^i(A) = 



cos 9k I — cos 9 m \ A { 1 + cos 0i, x 



cos 9 k I + cos 6*i, 



cos ir 



(C.3) 



(C.4) 



For large A the integrand [1 — U^ k (A)] forces k to stay away from the thrust axis 



by an small angle of order A 1 . Changing integration variables we have 



<D^_ / dA > ' ^ ' 1 



it^(A) 



dx 



x 



A' 



X 



A' 



1 -X 



cot 4 - X 



(C.5) 



At small A this gives the result (5.5), while at large A we obtain (5.6). 
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